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This review article aims at presenting the theory of inflation. We first describe 
the background spacetime behavior during the slow-roll phase and analyze 
how inflation ends and the Universe reheats. Then, we present the theory 
of cosmological perturbations with special emphasis on their behavior dur- 
ing inflation. In particular, we discuss the quantum-mechanical nature of the 
fluctuations and show how the uncertainty principle fixes the amplitude of 
the perturbations. In a next step, we calculate the inflationary power spectra 
in the slow-roll approximation and compare these theoretical predictions to 
\^5 ' the recent high accuracy measurements of the Cosmic Microwave Background 

f^ ' radiation (CMBR) anisotropy We show how these data already constrain the 

\l . underlying inflationary high energy physics. Finally, we conclude with some 

speculations about the trans-Planckian problem, arguing that this issue could 
allow us to open a window on physical phenomena which have never been 
probed so far. 
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1 Introduction 

Inflation is the most promising theory of the early Universe. It was invented 
H I by A. Guth pP at the beginning of the 80's in order to solve the puzzles of the 

hot Big Bang theory. A very interesting aspect of the inflationary theory is 
that it allows us to build a bridge between cosmology and high energy physics. 
This is particularly valuable in view of the fact that it is difficult to probe 
physics beyond the standard model of particular physics. 

However, the details of the underlying particle physics model are encoded 
into the fine structure of the cosmological observables. This is why, after the 
invention of the inflationary scenario and during quite a long time, it was 
in fact only possible to check the consistency of the inflationary predictions. 
The situation has now changed drastically with the recent releases of very high 
accuracy cosmological data. One can now take advantage of the full predictive 
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power of inflation with the hope to learn about physics in a regime which has 
never been reached before. 

The goal of this review article is to give a general presentation of the 
inflationary scenario. In particular, we will emphasize how the origin of the 
inhomogeneities present in our Universe is explained in the framework of in- 
flation. We will see that it is based on an elegant interplay between general 
relativity and quantum theory. Then, we will present the corresponding pre- 
dictions made by inflation and will study how the currently available data can 
already put some constraints on the underlying particle physics models. 

This article is organized as follows. In the next section, we describe the 
evolution of the inflationary background, the slow-roll phase and the reheat- 
ing. Then, wc present the theory of cosmological perturbations of quantum- 
mechanical origin and compare its predictions to the available data. Finally, we 
conclude this article with some speculations concerning the trans-Planckian 
problem of inflation, demonstrating that future astrophysical observations will 
maybe open a new window on high energy physics. 



2 The Inflationary Universe 
2.1 Basic Equations 

The cosmological principle implies that the Universe is, on large scales, homo- 
geneous and isotropic. As a consequence, the metric tensor which describes 
the geometry of the Universe is of the Friedman-Lemaitre- Robertson- Walker 
(FLRW) form, namely 

ds^ = -c^dt^ + a^{t)jlfdx'dx^ = 0^(77) l-dr]^ + j^fdx'dx^] , (1) 

(3) 
where 7^ ■ is the metric of the three-dimensional spacelike sections. The three- 
dimensional sections have a constant scalar curvature. The variable t is the 
cosmic time while 77 is the conformal time. They are linked by the relation 
cdt = adi]. In this article, we will work with dimensionless coordinates x^ and, 
as a consequence, the scale factor 0(77) will have the dimension of a length. 

The matter is assumed to be a collection of A^ perfect fluids and therefore 
its stress-energy tensor is given by the following expression 

TV 

T^,„ = ^ T^'J = (Pt + P-vhiiUi, 4- PT,gfi„ , (2) 

where p,^. is the (total) energy density and p^ the (total) pressure. These two 
quantities are linked by the equation of state, p^ — u {p^) [in general, there 
is an equation of state per fluid considered, i.e. pi = uji{pi)]. The vector u^ 
is the four velocity and satisfy the relation u^u^ = — 1. This means that one 
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has u^ — (1/a, 0) and u^ = (— a,0). The fact that the stress-energy tensor is 
conserved, ^"^T^^ = 0, amounts to 

p;+3-(p,+p,) = 0. (3) 

^ a 

This expression is obtained from the ^ = component. The component /z = i 
does not lead to an interesting equation for the background. If one assumes 
that each fluid is separately conserved then the above equation is valid for 
each species. 

We will assume that gravity is correctly described by the theory of General 
Relativity even in the very early Universe. This means that the equations 
which link the geometrical part to the matter part are nothing but the Einstein 
equations 

Rtiu - -zRgiJ.u = nTf_,^ , (4) 

where k = 8ttG/c'^ — 87r/m^ . These equations, in the case of a FLRW Uni- 
verse, are differential equations determining the time evolution of the scale 
factor and read 
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where a prime denotes a derivative with respect to conformal time. In the 
following, we will use the definition Ti. = a' /a. The parameter fc = 0, ±1 
represents the curvature of the spacelike sections. If, in addition, the equation 
of state of the perfect fluids are provided, then we have a closed system of 
differential equations and, therefore, the evolution of the corresponding model 
of the Universe is completely specified. 

2.2 The Inflationary Hypothesis 

By definition, inflation is a phase of accelerated expansion, i.e. the scale factor 
satisfies |2] 

It is interesting to postulate that such a phase took place in the very early 
Universe because, in this case, one can explain many different seemingly para- 
doxical facts like, for instance, the horizon problem or the flatness problem. 
Because of the latter, from now on, we will put fc = in the Einstein equa- 
tions. More precisely, one can show the inflationary scenario is satisfactory 
if the number of e-folds, i.e. the logarithm of the scale factor at the end of 
inflation to the scale factor at the beginning of inflation is greater than 60 [2] , 

iV^ > 60 . (7) 
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More detailed arguments about the advantages of inflation can be found in 
Ref. |3] but, at this point, it is important to notice the foUowing three facts. 
Firstly, inflation is convincing because, by means of a single concept or hy- 
pothesis, one can solve many different problems. In this sense, inflation is an 
economical assumption. Secondly, as we will show below, inflation is falsifiable 
since it makes definite predictions that we will describe. Therefore, there is the 
hope either to confirm or to exclude this hypothesis and, in any case, there is 
the certainty to learn something about the early Universe. Thirdly, inflation 
is defined by the condition Q but this does not prejudge the physical nature 
of the matter responsible for the acceleration of the Universe. The only thing 
one can say is obtained by expressing the acceleration of the scale factor in 
cosmic time. Using the Einstein equations, one gets 

- = -^(p + 3p), (8) 

where a dot denotes a derivative with respect to the cosmic time. Therefore, 
the fluid responsible for inflation must be such that 

P<-f, (9) 

i.e. must have a negative pressure. As a consequence, this fluid cannot be a 
standard fluid, like a gas for instance, but must be somehow "exotic" . However, 
this does not come as a surprise since inflation is supposed to take place 
at very high energies. At those energies, the natural description of matter 
is (quantum) fleld theory. As we are now going to demonstrate, it is quite 
interesting to remark that the most simple example of a field theory can do 
the job very well and "produce" the negative pressure which is necessary to 
inflation. We now discuss this point in more details. 

2.3 Implementing the Inflationary Hypothesis 

The most simple implementation of the inflationary scenario is to assume that 
matter is described by a scalar field (p{i]) 0E]. This case is nothing but a 
particular example of a perfect fluid. The corresponding action reads 



^g^'-d^ipd^ip + Vi^) 



S = - d^xV^ 

Then, the stress-energy tensor, which is defined by 

5S 



(10) 



2 SS 

T,,^~^=^-j^, (11) 



can be written as 



lg-f^d^^df!V + V{^) 



(12) 
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From this expression, it is clear that the scalar field is indeed a perfect fluid. 
The energy density and the pressure are defined by T°o = — p, T'^j = pS^j and 
we obtain 

The conservation equation can be obtained either by re-deriving it from the 
very beginning or just by inserting the previous expressions of the energy 
density and pressure into Eq. ||21l. Assuming ip' ^ 0, this reproduces the 
Klein- Gordon equation written in a FLRW background, namely 

^" + 24' + a^^=0. (14) 

a dip 

The other equation of conservation expresses the fact that the scalar field 
is homogeneous and, therefore, does not bring any new information. Finally, 
a comment is in order about the equation of state. In general, there is no 
simple link between p and p except when the kinetic energy dominates the 
potential energy, where oj = p/p ~ 1, i.e. the case of stiff matter or, on the 
contrary, when the potential energy dominates the kinetic energy for which 
one obtains w ~ — 1. This last case is of course very interesting since this leads 
to an inflationary solution. We have thus identified the condition under which 
inflation can occur: the potential energy must dominate the kinetic energy, 

We now turn to a systematic study of this regime. 



2.4 Slow-roll Inflation 

Since the kinetic energy to potential energy ratio and the scalar field accel- 
eration to the scalar field velocity ratio are small, this suggests to view these 
quantities as parameters in which a systematic expansion is performed. The 
slow roll regime is controlled by the three (at leading order) slow-roll param- 
eters defined by: 

The slow-roll conditions are satisfied if e and 5 are much smaller than one and 
if ^ = 0(e^, i5^, eS). It is also convenient to re-express the slow-roll parameters 
in terms of the inflaton potential. Using the equations of motion in the slow- 
roll approximation, one can show that 
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where, here, a prime denotes a derivative with respect to the scalar field. 

The equations of motion, that is to say the Friedman equation and the 
Klein-Gordon equation can be re-written exactly as 

H- = J^ ^ = ^^ (19) 

from which one deduces that, if the slow- roll conditions are satisfied 

H».|vw + oW. 'i--~^o(sy ,20, 

These equations are of course easier to analyze and solve than the original 
ones. 

Let us now analyze a concrete example. We choose the following class of 
potentials 

3A„ 4 f If 



where n is a free parameter and A„ the coupling constant. The factors that 
show up into the definition of the potential have been chosen for future con- 
venience. Let us first try to see under which conditions the slow-roll approxi- 
mation is valid. We adopt the criterion e < 1 (it is in fact e ^ 1 and , strictly 
speaking, e < 1 only corresponds to the condition necessary in order to have 
an accelerated expansion). This amounts to 

Tl 

4v7r 

Of course, this constraint applies in particular to the initial value of the field. 
We already conclude, that for this class of models, the values of the field must 
be at least a few Planck mass. Let us be more precise and evaluate the total 
number of e-folds during slow-roll inflation. It is given by the formula 

N^ ^ iJ^ c. -« r" d^F(^) (^) ' , (23) 



Vii 



dip 



from which one gets 



n \ m„, / 4 



Let A^min the minimum number of e-folds required in order to solve the prob- 
lems of the hot big-bang model (we have seen before that N^^^ — 60) then 
one has 

'^ini > TOp, J^ (iV„,i„ + ^) . (25) 
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For n = 2, this gives ip-mi — S.lm^-^ and for n — A, one obtains ip-mi — 4.4TOpj. 
However, it is often argued that "natural" initial conditions (see the last article 



in Refs. |21) are such that V{(pini) 



which amounts to 



Vh 



i/n 



A-i/'^ 



> TOp 



(26) 



because, as we will discuss later one, the COsmic Background Explorer 
(COBE) normalization implies that the coupling constant is small. In this 
case, the number of e-folds is a large number, much larger that the minimum 
required. Of course, the fact that the value of the scalar field must be larger or 
of the order of the Planck mass has led to many discussions about the model 
building problem. In this review, we do not address this question. Details 
about this issue can be found in Refs. |2|S|. 

Let us now solve the equations of motions in the slow-roll approximation. 
For the scalar field straightforward calculations lead to (n 7^ 4) 



ip{t) = ipi, 



n(4 — n) y/Xr, 






(4-")/2 



'71pi(i-iini) 



2/(4-«) 



(27) 



The last expression can also be expressed in terms of tend, the time at which 
slow-roll infiation stops. One obtains 



ip{t) = (fiu 



t-ii, 



'^cnd ''i: 



V V'ini 



(4-n)/2' 



2/(4-n) 



(28) 



The advantage of the above equation is to show that for t < t^nd the argument 
between braces always remains positive and hence the whole expression well- 
defined. A negative argument would simply signal the break-down of the slow- 
roll approximation and, in this case, the above formula cannot be used. 

Let us now turn to the scale factor. Integrating the first of Eqs. H2()(l leads 
to 

a{t) = flini exp \ -3:^ yiit) ~ (pL] \ . (29) 



nmt: 



From this expression, one can also calculate the evolution of the scalar field 
in terms of the number of e-folds N which is the natural time variable during 
inflation. One gets 



"'"'-"'"l/lSf -T."- 



(30) 



from which one obtains the formula giving the evolution of the Hubble pa- 
rameter during inflation 



HiN)^m^,VX„ 
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This equation is valid until A^ = N^ and, in this regime, the above for- 
mula is always well-defined. Indeed, Eq. (|^ becomes meaningless at A^max = 
47r((/3ini/mp|)^/n but -/Vmax > -^T- The above equation has interesting conse- 
quences for our understanding of inflation. It shows that the Hubble parameter 
can evolve and change significantly during the slow-roll phase. Later on, we 
will see that a quantity which plays an important role is the value of the 
Hubble parameter when the scales of astrophysical interest today crossed out 
the horizon during infiation. This happens 60 e-folds before the end of infia- 
tion. This scale is constrained by the observations on the Cosmic Microwave 
Background Radiation (CMBR) anisotropics to be H/m^^ < 10~^. However, 
this does not mean that the Hubble parameter has not been larger before, 
especially if the total number of total e-folds is large, as it is the case for the 
initial conditions discussed in Eq. I|26() . For instance, if we have a massive po- 
tential, n = 2, and ipini — lOOrripj then inflation starts with an initial Hubble 
parameter of Hi^i ~ 10~^mpj but ends at ilond = Wpi%A^ [?^/ (4\/7r)]' ~ 
0.28 X 10-^mp, after N ~ 62000 e-folds, where we have used A2 ~ 10"^° 
(corresponding to a mass m ~ 10~^mpj). Therefore, in summary, it will be 
important to keep in mind that the observations give indications about the 
scale of inflation when the relevant scales crossed out the horizon during in- 
flation but cannot, a priori, put constrains on the Hubble parameter in the 
earliest phases of evolution. This situation is summarized in Fig. ^ 

When the field reaches the value tpcnd, slow-roll inflation stops and the 
system enters a new regime that we now briefly describe. 

2.5 Reheating 

When the scalar field reaches the point where the slow-roll parameter e ~ 1, 
for which ip = ifend, inflation stops and the flcld starts oscillating around its 
minimum OC]. In this regime, the system is governed by two time scales: 
the Hubble time H~^ and the period of the oscillations around the minimum 
{V")~^ (here, a prime means derivative with respect to the fleld). The im- 
portant point is that these two scales are very different. The frequency of the 
oscillations is much larger than the Hubble rate, Wosci — V" ^ H. The scalar 
field obeys the Klein-Gordon equation which can be put under the following 
form 

-3Hip^ ^~6H{p-V), (32) 



dt 

where the relation ip"^ = 2{p — V) has been used. This equation can be time 
averaged and one gets 



-{6H{p-V)) o:i-6H{p-V) , (33) 



where we have used the fact that the Hubble rate does not change during one 
period of the oscillations. The right hand side of the equation above can be 
evaluated as 
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Fig. 1. Evolution of the various scales discussed in the text during inflation and 
the subsequent radiation and matter dominated epochs. In particular, it is apparent 
that the CMBR measurements only probe the inflationary model when the modes 
of astrophysical interest today crossed out the horizon during inflation. The small 
window shows a typical inflationary potential, see Eq. (I2H . 



\fp^V{^Aip 



Akp 



^P~V{^) 



{P~V)^-J {p-V)dt = 

(34) 

where (pm is the value of the scalar field at the maximum of its oscilla- 
tions. In order to obtain the previous relation, one has also utilized that 
di — A^pI \J'l{p — V). Then, one uses that over one period, p ~ V{<Pvn) = Vm 
is a constant and one obtains 



(p-V) -IP, 
where the number 7 is defined by 
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the last result being valid for potentials of the form V{ip) oc (/?". Let us now 
turn to the left hand side of Eq. (|33|l . The term {dp /At) can be written as 
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Ap/T . It can be expressed as a finite difference expression and one can rewrite 
it as p. Tliis is valid for time intervals much larger than the period of the 
oscillations. Therefore, the equation governing the evolution of the energy 
density of the field H33|) can be rewritten as 

p = ^i/p^poca-6"/("+2). (37) 



Then, the scale factor is given by a{t) ex t("+2)/(3"). For the massive case, 
n = 2, the energy density evolves as in a matter-dominated epoch. This can 
be easily understood since, in this case, the Klein-Gordon equation is exactly 
the equation of an harmonic oscillator. In this situation, it is known that 
((/3^/2) = {V{ip)) which implies that the pressure vanishes. 

So far, we have not taken into account the effect of particles creation. 
Phenomenologically, it can be described by adding a term Fip in the Klein- 
Gordon equation which now reads 

^ + 3i/^ + r^ + ^^o ^^ = -^(3i7 + r)p. (38) 

dip at n + 2 

If we assume that the particles produced are very light in comparison with the 
mass of the inflaton, these particles will be very relativistic. This means that 
the equation of conservation of radiation should also be modified according to 

^ = -mp, + rp, (39) 

so that the total energy is still conserved. Eq. (|S5|l can be easily integrated 
and the solution reads 



-6n/(n+2) 

pit) = Posci ( —^— 1 exp 



2n 



■r{t-t. 



(40) 



where t = iosci is the time at which the oscillations start (i.e. the time at 
which the slow-roll period ends) and posci is the value of the scalar field energy 
density at that time. The effect of the term Fip is to multiply the result H37|) 
by a decreasing exponential factor. Equipped with this solution, we can solve 
Eq. (|39|l and determine the evolution of p,-. If the scalar field energy density 
dominates the radiation, as it is the case at the beginning of the reheating 
period, the solution reads (using the fact that the scale factor is known in this 
regime, see above) 



/ \ -4 / , r, \ (n-|-8)/(3n) / „ 

I a \ I n + 2 \ I 2n 

p,-{t) = TiosciPosci TT-r^ exp — — -riosci 

V aosci / V 2?ii iosci / V "- + 2 

fn + 8 2n \ fn + 8 2n 



(41) 



where the function 7(0, x) is the incomplete gamma function defined by 
7(a, x) = J„ e~*i"^^dt. In the above formula, we have assumed that, at the 
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end of the slow-roll epoch t = iosci: Pi- — 0. For small values of the argument 
X, the incomplete gamma function reduces to ~ x" /a. We define t^^^ = F^^ 
and then we have x = 2nrt/{n -I- 2) = 2n/{n + 2){t/t^^) and for times 
t > iosci ^ r~^, the argument of the incomplete gamma function is small. In 
this limit, one obtains 



Pr[t) ^ 1 Poscit, 



^osci'-osci 



{n + 8)t 



1- , 

tr. 



^ \ -{n+a)/{3ny 



(42) 



We see that p,. starts to increase, reaches a maximum and then decreases. 
When t approaches t^^^ , the previous approximation breaks down since the 
argument of the incomplete gamma function is no longer small. However, for 
order of magnitude estimates, we can try to push this approximation. At i ~ 
tj^jj , we have pr — r poscitosd'^n / [{n+8)t-^-^] since the second term in the square 
bracket in Eq. (|42|l is negligible. After thermalization, the energy density of 
radiation takes the form p^ — g^Tr^T^/SO. Using the fact that posci — H^^fin^^ 
(nothing but the Friedman equation) and that iosci — ^hif ; '^^^ can deduce 
the reheating temperature 

^- - -^9. [^^^) irm^.) ■ (43) 

Then, from this temperature, the universe evolves in a standard radiation 
dominated era. The remarkable feature of the previous equation is that it 
does not depend on the scale of inflation H^^^ but only on the decay rate F of 
the inflaton. This means that whatever the scale of inflation is, the radiation 
dominated era always starts at the same energy (at fixed decay rate) and that 
the duration of the period of coherent oscillations can change quite a lot. The 
number of e-foldings during this epoch can be evaluated as [since during this 
epoch, the scale factor scales as ex t("+2)/(3")] 

The previous considerations are valid if the life time of the inflaton is bigger 
than the age of the universe at the end of inflation. Otherwise, there is no 
period of coherent oscillations. In this case, the vacuum energy H?^fm^^ is 
directly converted into radiation and the reheating temperature is 

OQl/4 

T,,^^g:'/\H,^,m,,f^\ (45) 

Finally, let us recall that the calculations above assume that the physical 
quantities are time averaged and therefore that the time scales considered 
are larger than the period of the oscillations. In Fig. El where the evolution 
of the fleld versus the number of e-folds is displayed, we have integrated the 
equations of motion numerically. This plot conflrms our analytical estimates: 
inflation consists of a phase of slow-roll followed by a phase of oscillations. 
This concludes our study of the inflationary background. 
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Fig. 2. Evolution of the scalar field during slow-roll inflation and the reheating phase 
(where the field oscillates) obtained by numerical integration of the equations of 
motion. The potential is of the type of Eq. l|^ with n = 2 and A2 = (87r/6) x 10~^°. 
The initial conditions are such that ipini — Snip^ leading to A'^^ ~ 60 as confirmed 
by the plot. 



3 Cosmological Perturbations 
3.1 General Framework 

It is an observational fact that the universe is not isotropic and homogeneous. 
Therefore, if one wants to have an accurate description, it is clearly mandatory 
to go beyond the FLRW model. On the other hand, it is also an experimental 
fact that, in the early Universe, the deviations from the isotropy and from 
the homogeneity were small (e.g. from the COBE measurement, ST/T ~ 
10~^). This suggests a perturbative treatment. Therefore, the following metric 
tensor |S] gives a refined description of our Universe 



It^i^iV, X) = [gf,,y{v) + <^hf,^{T], X) + 



e%^(?7,a;) + 



■]dx^'dx'' 



(46) 



where gfj.i^irj) is the standard FLRW metric introduced previously and repre- 
sents the "background" (the parameter e in the above equation should not be 
confused with the first slow-roll parameter; they have nothing to do with each 
other). The perturbed metric depends on x and this is the signature of the fact 
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that we now go beyond the cosmological principle. In order to be consistent, 
the same expansion must be performed for the quantities describing matter. 
For example, if there is a background scalar field ip{r]), a refined description 
of the scalar field can be expressed as 

(p{T], x) = Lp{rj) + e&if{r], x) + e'^S^'^'>ip{T], x) + ■ ■ ■ . (47) 

The main goal of the theory of cosmological perturbations is to determine the 
evolution of the perturbed quantities /i^^ and Stp and, then, to use them in 
order to calculate observable quantities. To find the behavior of the perturbed 
quantities, one needs some equations of motion. Naturally, these equations 
are taken to be the perturbed Einstein equations written order by order (we 
assume that gravity is described by General Relativity). Therefore, we expand 
the Einstein tensor and the stress-energy tensor according to 

G,.. = G|,"j + 6G|,ij + e^Glfj -I- . . . , T^, = T^^J + eT^^-J + e^T^^J + • • • , (48) 

and then identify the terms of same order to obtain 

G^v — ^^/ii/ ' G^iy = '^^/jiy ' G^y = ''^T^^ , • • ■ • (49) 

In the present context, we will restrict ourselves to the linear order in the 
parameter e. 

Let us now try to describe the perturbed metric tensor in more details. 
For any symmetric two-rank tensor, there is a theorem [S] which states that 
hfj,,y{ri,x) can be decomposed as h^^{ri^x) — h\j_J + h]iJ + hl^J , where /i/ii; is 
constructed only from scalar functions, hliJ is constructed only from three- 
dimensional vectors with vanishing divergences and hliJ is obtained only from 
transverse and traceless three-dimensional tensors. These three types of per- 
turbations are the scalar, rotational and tensorial fluctuations respectively. 
Explicitly, the theorem implies that the unperturbed metric plus the per- 
turbed metric can be expressed as |S] 

ds^ = a^{ri){-{l + 2(j})ATf + 2{d,B - S'^dxMr; + [(1 - 2ij)5^.j + 2d,djE 

+djF, + d,Fj + h'-^^]dx'dx'} , (50) 

with Si and Fi being transverse vectors, i.e. d^Si — d^Fi — and hl^ being 
a transverse and traceless tensor, i.e. S^^hij — 0, d^hij = 0. We see that the 
scalar part of the metric depends on four unknown functions: 0, B, "0 and 
E. The vector part depends on two vectors with vanishing divergence, i.e. Si 
and Fi and, finally, the tensor part depends on one transverse and traceless 
tensor, namely /i^ . At linear order, each type of perturbations decouple and, 
as a consequence, can be treated separately. 

In the specific case of inflation, one can show that vector perturbations 
cannot be produced |S]. Hence, in the following, we will consider that only 
scalar and tensor perturbations are present. 



14 Jerome Martin 

At this point, one should discuss a weU-known problem of the theory of 
cosmological perturbations: the gauge issue. A complete study of this ques- 
tion can be found in Refs. [Sl llUlITT] but, roughly speaking, it consists in the 
following. There exist solutions to the perturbed Einstein equations which 
are coordinates dependent, i.e. which can be removed by performing an in- 
finitesimal change of coordinates. These solutions are fictitious and should not 
be considered as physical. The following analogy may help to understand the 
problem jEl ■ Let us consider the four-dimensional FLRW manifold denoted V4 
in what follows. It can be embedded into a higher dimensional manifold, more 
precisely into the five-dimensional Minkowski spacetime Ef 4 whose metric is 
Vab where the indexes A and B runs from to 4. A point in Ei 4 is located 
by its coordinates z^. An embedding is a map from V4 to E^ ^■. z^ = z^{x'^). 
For a spatially fiat FLRW spacetime endowed with Cartesian coordinates, the 
embedding explicitly reads: 

z°{r,,x,y,z)^^aiTj)ix^+y^ + z^ + l) + ^J ^TT^^lr, (51) 

zi(,7, X, y, z) = \a{^){x^ + y^ + ,^ - i) + ^ T ^^dr , (52) 

2 2 J a'[T) 

z'^{ri,x,y,z) — ax, z^{r],x,y,z) — ay, z {■q,x,y,z) — az . (53) 

Therefore, the FLRW manifold can be viewed as a surface into the higher 
dimensional spacetime E\^. The metric of this surface can be calculated by 
means of the well-known formula 

Qf^i, (77, x) = riABdf,z^dyZ^ , (54) 

and we can indeed check that this reproduces the metric of a spatially fiat 
FLRW universe. Let us now try to "deform" this manifold since this is what we 
have in mind when we consider small perturbations around the background. 
In the present context, a deformation consists of the following. If we consider 
a point M in the manifold V4 located by its coordinates z^{x^) in E^ ^, 
deforming the manifold means slightly displacing the point M in E^ ^. This 
means that the coordinates of this point are no longer z^ but z^ + ev"^{x^) 
where e is a small parameter. The vector v {x^) characterizes the deformation. 
As a consequence, the new metric of V4 calculated by means of Eq. H54I) reads 

3m"^ = 7^1^ + '2er]ABdf,z'^d^v^ . (55) 

However, all the vectors v^{x'^) do not represent a genuine deformation. In- 
deed, if the following relation is satisfied 

z'^ix'') + ev^ixf") = z'^ix'' + e^) , (56) 

then, clearly, the displacement is within V4 and does not correspond to a defor- 
mation. This is merely a change of coordinates that should not be considered 



Inflationary Cosmological Perturbations of Quantum-Mechanical Origin 15 

as a physical deformation of V4. This gauge problem consists of identifying the 
spurious modes and in removing them from the theory. To conclude this di- 
gression, it should be emphasized that the link between the previous approach 
and the theory of cosmological perturbations has never been worked out in 
details. Therefore, an important warning is that it may well turn out that the 
analogy used above cannot be applied completely to the theory studied here. 
Having realized that there are non physical modes, the problem is now 
to find a method to get rid of them. Following Bardeen's seminal paper, an 
efficient way is to work with combinations of the metric tensor components 
which are invariant under a general change of coordinates (a "gauge" trans- 
formation) and, hence, which cannot contain a spurious mode. For scalar 
perturbations, the two following combinations ^01 



^{ri,x) = (j)+ - 



a{B-E') 



W{ri,x) = ^--{B-E'), (57) 



a 



are gauge invariant. They are called the Bardeen potentials. In what follows, 
we will see that, in the simple case where matter is described by a scalar field, 
one has in fact <I> — 'F. This means that we have reduced the study of the 
scalar perturbations to the study of a single quantity: the Bardeen potential 
^i],x). 

The case of gravitational waves remains to be treated. In fact, it is easy to 
realize that the gravitational waves are gauge-invariant by definition because 
one cannot construct an infinitesimal change of coordinates with a tensor. 
Therefore, one can safely work with the tensor h^- {rj, x) introduced before. 

We have identified the gauge invariant variables that describe the gravi- 
tational sector. Our next move is to do the same but for the matter sector. 
This can be done in general |ll)[lll| but, since we have inflation in mind, we 
just consider the case of a scalar field. Then, one can show that 

6ip^'^'\T],x) = S'fi + 'fi'{B^E'), (58) 

is the gauge-invariant perturbed scalar field. 

Finally, we need a last ingredient. Since the spacelike sections are fiat and 
since we study the linear theory, it is very convenient to work in the Fourier 
space. Indeed, because of the above properties, each Fourier mode evolve 
independently (the mode coupling appearing at quadratic order only) and 
it is sufficient to follow their time evolution. Therefore, we Fourier transform 
the Bardeen potential and the gravitational waves according to 

HV, x) - j^^ I dfe<?(ry, fe)e''=- , (59) 

1 

(M 



C(^'^) - 7^ /dfc E PUkKiv,k)e^''-'- (60) 

^ ' s=+,x 



In the last equation, Pij{k) is the transverse and traceless polarization tensor 
satisfying the following properties: p^Ak^p^^" (fc) — 25"^ . The symbols "-I-" 
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and "x " denote the two possible states of polarization of the gravitational 
waves. Of course, we also Fourier transform the perturbed scalar field and 
work with Sip^^^'{ri, k). 

Having identified what the relevant degrees of freedom are, we now turn 
to the question of establishing their equation of motion. 

3.2 Equations of Motion 

Since the Einstein equations are obviously gauge-invariant "by definition" , it 
is clear that it is possible to express them in terms of gauge invariant quan- 
tities only. We start with density perturbations. Lengthy but straightforward 
calculations lead to (for a fixed Fourier mode fc) 



-3nin'P + ^') -fc^if = - 



Ti'P + W' = -ip'Sip'-'^''^ 



Wf^- 



w = o, 



Sip 



(gi) 



2dV 

a - — dip 
dip 



(gi) 



{2H' + H^)<P + H<P' + W" + 21-1^' - \k^{<P - •f ) = ^ \-{lp')^<P 



+ip' S(p 



o(gi) 



dip 



,(61) 
(62) 

(63) 



As announced, the two Bardeen potentials <P and ^ are equal. This is true as 
long as there is no anisotropic stress. Despite the apparent complexity of this 
system of equations, straightforward manipulations show that everything can 
be reduced to the study of a single equation which reads 



<P" + 2\H 



p' 



$' 



2\n' -n 



p' 



<? = o. 



(64) 



This equation is valid provided p' ^ 0. In this case, for which the scalar field 
plays the role of a cosmological constant, we have no density perturbations at 
all, <? = 0. This does not mean that the perturbed scalar field cannot fluctuate 
in de Sitter spacetime (as a matter of fact, it does) but that, in this case, these 
fluctuations do not couple to the fluctuations of the metric. Eq. (|64ll can be 
transformed in order to permit a more transparent physical interpretation. If 
we consider the variables u and Q defined by 



u = ^ 

3 n 



p 



1/2 



V3. 



ap}' 



3 1 



2a^- 



(65) 



then Eq. H64|l can be expressed as: 



(66) 
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The above equation can be viewed either as the equation of a parametric 
oscillator, with a time-dependent frequency given by uj^{k,ri) = fc^ — 0"/9, 
or as a Schrodinger equation with the potential 9" /9. This effective potential 
contains derivatives of the scale factor up to the fourth order. The typical 
behavior of the solutions can be easily found. For modes k^ 3> 9" /6, the 
variable u oscillate, u ex e^^"^ , while for modes fc^ <C 9"/ 9 the solution of 
Eq. H66|l may be expanded in powers of fc^ . At leading order we obtain 

u(r;, k) = A^{k)9{r^) J' -^ + A2{k)9{^) . (67) 

Since — > oo for a — > in general, Ai is the arbitrary constant in front of 
the regular (growing) mode and A2 a constant associated with the singular 
(decaying) mode. We will see in the following that the variable u is in fact 
not the most interesting for density perturbations. Quantum-mechanical con- 
siderations, among others, will lead us to work with another variable. From 
the above solution for u(?7, fe), we easily deduce the Bardeen potential in the 
superhorizon regime. One obtains 

<?(r;, k) ^ ~A,{k)^f a'j{r)dr + ^J ■ (68) 

For example, for a power-law behavior of the scale factor, i.e. a ex JTyj^^'^ with 
P < —2 in order to have inflation, the 'growing' mode turns out to be constant 
in time, namely 

<P{v,k)c,^j±^A,ik). (69) 

For the de Sitter case, w = — 1 and we recover the fact that there no density 
perturbations at all. 

The equation of motion H64(l has a first integral for modes that are much 
larger than the Hubble scale, i.e. k/a ^ H. Following Ref. |H] we define 

which was introduced by Lyth |13| originally. Essentially, the quantity ^ is 
the perturbation of the intrinsic curvature in the comoving gauge |13| and 
is written —'R. in that reference. The equation of motion for the Bardeen 
potential can be re-written as an expression for the first derivative of the 
quantity C- One obtains ^J 

Of course, to derive this equation we have not assumed that the equation of 
state parameter is a constant. Thus, C is constant in time for superhorizon 
modes kjH. <C 1 since then C = 0. 
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The importance of the quantity ^ is due to the fact that this is a pure 
geometrical quantity. Concretely, this means that the conservation law estab- 
lished above in the case of a scalar field is in fact valid for any type of matter 
(at least, provided that the so-called entropy perturbations do not play an im- 
portant role). Therefore, ( can be used as a "tracer" of density perturbations 
regardless of the type of matter responsible for those fluctuations. In particu- 
lar, it can be used to propagate the spectrum from the end of inflation (where 
the Universe is dominated by a scalar field) to the radiation dominated era 
(where the Universe is dominated by a relativistic fluid) without knowing the 
details of the reheating process. Let us now study how the calculation works 
in details. On superhorizon scales, the Bardeen potential is almost constant. 
If we neglect its time derivation then the equation of motion for (p and the 
definition of C, lead to 



$ 



3(l + ^.nf) ^ ^ ^x,Jei) 



3uj. 



< 



2H 



Sip"- 



(72) 



Using the equation of motion of the background, the above formula can be 
put under the following form 



C(^,fc) = 



3llj. 



'-n 



,5^fe')(r;,fe) 



f' 



(73) 



Now, let us assume that we want to know the Bardeen potential in an era 
dominated by a fluid with a given equation of state uj (concretely, we have 
in mind lo — 1/3 or lu = for the radiation or matter dominated epochs, 
respectively). Writing the constancy of C, i.e. Q^^: = C^^, one arrives at 



$Uv,k)^3n 



1 



5 + 3w 



V^')(?7,A;) 



(74) 



where we have used u^^j ~ — 1. This equation is very important since it links 
the primordial fluctuations of the scalar field to the fluctuations of the gravi- 
tational potential during the subsequent phases of evolution of the Universe. 
Let us now turn to gravitational waves. In order to obtain the equation 
of motion, we must compute the perturbed Ricci or Einstein tensors for the 
metric given in Eq. H50|l. One finds that SR'^o — 5R^i — 0. This result is con- 
sistent with the fact that the gravitational waves are transverse and traceless 
since only the trace and/or the derivative of the metric tensor can appear 
in these components. On the other hand, the component dR^j, i y^ j is non 
vanishing and the leads to 



SR' 



1 

2a2 


"/iW^; 


" a' 

+ — 


"/i(T)»; 



1 

2^ 



dkd'^h^^^'j = kST'j , (75) 



where 6T^j represents the anisotropic pressure part of the perturbed stress- 
energy tensor. For a perfect fluid (e.g. for a scalar field), it vanishes. However, 
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it is important to keep in mind that, a priori, the gravitational waves are not 
"independent" from the matter fluctuations and, hence, that they should be 
considered on the same footing as density perturbations. This is conceptually 
important because this means that it would be incorrect to argue that both 
types of perturbations should be treated differently, in particular with respect 
to the quantization of the cosmological perturbations. 

Then, for the rescaled Fourier amplitude defined by A*^(7y, k) = a{ri)h^{rj^ k), 
the equation 6R''j — can be re-written as J3] 



K) +(^fc --JM^=0. (76) 

This equation can be viewed either as the equation of a parametric oscillator, 
i.e. an oscillator with a time-dependent frequency, 0^^(77, k) = k'^ — a" /a or as 
a "time-independent" Schrodinger equation with a potential U,^(ri) — a" /a. 
Therefore, we obtain the same type of equation as for density perturbations. 
However, it is also interesting to notice that the effective potential for tensor 
perturbations involves the scale factor and its derivatives up to second order 
only. This difference is especially important during the reheating phase. 

As it was the case for density perturbations, it is clear that the solution 
to the equation of motion possesses two regimes. If the wave number k is 
such that fc^ » U^, then the mode function oscillates, i.e. ^^ ex e''^''. The 
interaction with the barrier (if any) corresponds to the time krjt ~ 1 since, in 
the inflationary phase, one has U^ ~ l/if (for a power-law scale factor; this 
is also the case for slow- roll inflation, see below). This time is also, roughly 
speaking, the time of Hubble radius exit. Indeed, the wavelength is given by 
A = 2TTa{r])/k and the Hubble scale is H^^ — a/H. The condition A = H~^ 
gives krjt ~ 1 since Ti. ~ I/77. The second regime is when the wave is below 
the potential, k^ <^Ur^. An approximate solution is 

dr 

where Bf(k) and -B|(fc) are two constants which are a priori free. The first 
term in Eq. H77|) is the growing mode whereas the second term is the decaying 
mode. This can be seen, for example, if we consider scale factors of the form 
a{rj) =4|r?|i+^. In this case, /x^ ~ Bf (fc)|?7|i+'5 -|-B|(fc)|r/|-^ and for /3 ~ -2, 
the first term goes to infinity while the conformal time goes to zero at the 
end of inflation. Therefore, this is indeed the growing mode. In terms of the 
amplitude h'^{ri, k) itself, one sees that the growing mode corresponds in fact 
to a constant and hence is conserved on large scales. Somehow, /i''(?7, k) plays 
for gravitational waves the same role as ( for density perturbations. 

Let us end this section with a comparison between density perturbations 
and gravitational waves. Using the results obtained before, the tensor to scalar 
amplitudes ratio today is given by 



<(,.»:) = b; (*:)<.(,) + B| (<:)<.(,) / ^^ , (77) 
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Therefore, if we assume that K^^ ~ ^i„fi which is the case if the perturbations 
are of quantum-mechanical origin, then we have h^/^u *^ 1; i-e. scalar fluctu- 
ations dominates over tensor fluctuations, because during inflation oj.^^^ ~ — 1. 
The previous considerations also illustrate the limitations of the classical 
approach. Without a theory of the initial conditions, i.e. without a prescription 
to choose the fc-dependent constants Ai^2{k) for density perturbations and 
Bi^2{k) for gravitational waves, we cannot really go further. This will be one 
of the main advantage of the quantum-mechanical version of the previous 
theory: a natural choice for Ai,2{k) and Bi^2{k). 

3.3 The Sachs- Wolfe Effect 

The production and the ampliflcation of small inhomogeneities in the early 
Universe described above has several observational consequences. In this re- 
view, we focus on one of them: the presence of small angular anisotropics in 
the temperature of the Cosmic Microwave Background Radiation (CMBR), 
at the level of ST/T ~ 10-^ detected for the first time by the CORE satel- 
lite in 1992 ^^. These anisotropics are of utmost importance for the theory 
of inflation because they allow us to check the predictions of this scenario 
and/or to constrain the physics of the early Universe. We now turn to a rapid 
discussion of this effect, a complete presentation being available in Ref. |16| . 

The Sachs- Wolfe effect jEl links the angular variations of the temperature 
on the celestial sphere to the presence of cosmological fluctuations in the 
early Universe. We have to calculate the change in the energy of the photons 
propagating from the last scattering surface to Earth. This energy is given by 
E — —^^vU^^k^, where k^ is the wave vector of the photon and u^ the velocity 
of the observer. Let us flrst investigate this relation for the background. 

Fundamental observers are observers who move with the cosmological flow. 
A trajectory is given by the set x^ = x^{s), where s is a afSne parameter along 
the line. The velocity along this curve is given by u^ = dx'^/ds and satisfles 
u'^Ufj, = — 1. For a fundamental observer, one has u' = by deflnition and the 
normalization of the four velocity implies that u'^ — (l/*^; 0) ^-nd u^ = (—a, 0). 
Let us now study the propagation of a photon. If fc^ = dx'^/dA is the wave 
vector of a photon, then the path followed by this photon is such that 



dfc^ 
"dA" 



ri^pk^kp = . (79) 



The solution of this equation is the trajectory of the photon: x^ = x'^{X). In 
addition, we have the constrain fc'^Tc^ = 0, expressing the fact that the photon 

follows a null geodesic. At zeroth order, this constraint gives Sijk^y = {kP^ . 
In the non perturbed universe, Eq. (|79|l possesses the solutions k^ — C^ /a?, 
where C° is a constant and /c* — —C^e^/a^^ where e' is a three vector such 
that de'/dA = 0. Taking the ratio of the wave vector components, we deduce 
that dx'/dTy — —e^ . Finally, integrating this relation, we find the equation of 
the trajectory in the unperturbed Universe 
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-e\v-Vu) + ^u 



(80) 



where (?7D,a;J^) are the coordinates at detection of the photon. It does not 
come as a surprise that the photons propagate along a straight hne. On the 
other hand, the energy is given by 



Eiv) 



0(77) 



(81) 



and, in fact, we just recover the well-known time evolution of the temperature 
(which, as expected, does not depend on space for the unperturbed Universe). 
Let us now turn to the Sachs- Wolfe effect itself. Essentially, this consists in 
computing the energy of the photons at first order. In a perturbed Universe, 
the most general observer possesses a velocity given by u'^ -I- Su^^, where u'^ 
denotes the velocity of a fundamental observer calculated above and where we 
assume that the components of 6u^ are small with respect to this fundamental 
velocity. The fact that the total velocity is normalized to —1 implies that 
6u'^ = —(j)/a. We also write Su^ as 6u^ = v"^ /a, from which we deduce that 
6ui = avi + adiB. The trajectory of the photons can also be expanded as 
a;'' -|- 5x^, where x'^ is the path of the photon in an unperturbed Universe 
determined before and 5x^^ are the small corrections around the background 
trajectory due to the presence of the fluctuations. In the same manner as 
we did for the four-velocity, we can expand the wave vector of the photon 
according to fc^ + 5k^^, where 6k^ = d {6x^^) /dA. At first order, the variation 
of energy can be expressed as 



5E = -h^^u^'k" - g^^Ju^k" - g^.yU^'Sk" , 



(82) 



which can be re- written as 6E = C'~'(l)/a + C'^e^{diB + Vi)/a -f aSk'^ . In this 
equation, the only unknown quantity is 5k^ and we now establish its expres- 
sion. Integrating the perturbed version of Eq. H79|) . one finds that 



<5fc" 



C" 



a^iv) Jr, 



dr 



2e'a. 



e'e^d.djB 



+h. 



(T) 



-2i}6, 



2d,djE 



(83) 



where rf-^ is the conformal time at emission. Let us stress again that the 
integration is performed along the unperturbed path of the photon. Putting 
everything together, we finally obtain 



§^ = ^^i+[0 + e^(a.s + ..)]:- 



dr 



2e'a, 



'e^d,diB 



iWjJ 



-2Vj(5, 



2didjE 



AT) 



(84) 
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The above expression depends on the coordinates of emission and detection of 
the photons. To go further, it is necessary to specify the conditions of emission, 
that is to say the characteristics of the last scattering surface. At zeroth order, 
the surface of last scattering has coordinates ry^, = 771SS, x^. = —e''{r]iss—'r]^)+x'^^ 
where rjiss is fixed and corresponds to the redshift ziss — 1100. The only de- 
pendence is now the vector e' and this corresponds to different directions on 
the celestial sphere. However, in presence of perturbations, emission occurs 
at different times and at different positions. In other words, the time of emis- 
sion is given by 77,^ = Tyiss -I- Sri{riiss,xl,). The quantity Sr]{r]iss,xl,) depends on 
our definition of the surface of emission. Let us assume that this surface is 
such that the density of photons, p^, is constant. Writing this condition at 
first order gives Sp^{r]iss, xl,) + p'^{i]iss)Sri{riiBs, x^,) = 0. Using the conservation 
equation which implies that p' — —AHp^, we arrive at 

oviViss.x ) = —— — . (85 

"" 47^(r]ks) Pjiviss) 

Therefore, the term a{ri^) in Eq. 1841) should be written as 

"'{Ve) = a(77iss) +'H{itiss)St]{t]iss,xI.) = 0(77133) + -(5^,(77133,2;^) , (86) 

where S-y = Spj/ pj is the density contrast. In the same manner, if we say 
that detection takes place on a surface such that the baryons energy density 
is constant, the factor a^^(7/n) should be written as 0,^^(1]-^) — a^^(7/o)[l — 
{l/3)Sh{'rio,x'^ )] (the factor 1/3 comes from the equation of conservation but 
now written for a fluid whose equation of state vanishes). Finally, Eq. H84() 
takes the form 

^-^{1 + ^-5,(77133,0 -i^b(%,<)+[0 + e'(a,S + t,)]: 



dr 

»7laB 



■' - 2e'd^cf) - e'e^d.djB + -^jS.j + d^d,E 



2% 



(87) 



Having established this important relation, we must now show that this ex- 
pression is gauge- invariant. For this purpose, it is sufficient to express the ra- 
tio of the energies at emission and detection only in terms of gauge-invariant 
quantities. We have already described the gauge-invariant variables for the 
gravity sector. For the variables describing matter, we only need the gauge- 
invariant density contrast 5g = 5 + p' / p{B — E'). Finally, one has to decom- 
pose the three- velocity as Vi = diV and the gauge- invariant velocity can be 
expressed as v'^^^' = v + E' . Let us also notice that the spatial derivatives 
can be expressed in terms of time derivatives. Indeed, along a trajectory, one 
has df/drj = djjf — e^dif from which we find e^dif = 0,,/ — df/drj. Then, 
straightforward calculations show that 



E^ aim) { 4--^^^'-^'-^-^ 3^ 
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We have thus proved the gauge invariance of the ratio E^ /E^ TT^ . 

The Sachs- Wolfe effect is frequency independent. This means that the 
shape of the black body is preserved at the perturbed level and this is why 
a perturbed temperature is still a meaningful concept. If we define ST/T = 
[STj^ — Tol/^D with Tq = Tjja(?7iss)/a(?7o), we arrive at the final form of the 
Sachs- Wolfe effect, namely 



with, 



ST\ 



T \T I \T I \T I 



(D) 



\t) ~ 4^"^ 



((5^)g(?7k.,xM +<Z'(7yk.,a;M - e'^^v^^\^J^,,,xl) 



fVo 

/ dr {<P' + If') 



Several comments are in order here. Firstly, we have discarded the terms 
#(r?o,x' ) and 5^^ (tjqjX"^ )/3 since they do not depend on the vector e'. Sec- 
ondly, the first term [ST/T]^^") has its e* dependence fixed. This is just the 
dipole term due to our motion with respect to the frame of the CMBR. 
Thirdly, the other terms are genuine fluctuations of primordial origin. As 
already mentioned, they have been discovered in 1992 by the COBE satellite. 
Finally, let us conclude this section by establishing the expression of the 
Sachs-Wolfe effect due to density perturbations on large scales. On these 
scales, the Doppler term is negligible. The integrated Sachs- Wolfe effect is 
also negligible because, on superhorizon scales, the Bardeen potential is ap- 
proximatively constant, hence its derivative vanishes (see before). Therefore, 
only the first two terms remain. One can show that they combine such that 

— j ^-<?(ryi,„<). (91) 

This equation permits to compute the angular power spectrum in the COBE 
regime, i.e. for large angular scales. 
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4 Quantization of Cosmological Perturbations 

We start this section with a discussion of the quantization of a free scalar 
field. This constitutes the prototype of methods used in the sequel for the 
cosmological perturbations. 

4.1 Quantization of a Free Scalar Field 

We consider the question of quantizing a (massless) scalar field in curved 
space-time. The starting point is the following action 

S ^ --^ j A^x^—gg^'']^d^,$d,$ , (92) 

which, in a FLRW Universe, reads 

S=Y I ^^xa^iv) i<P'^ " S'^d^^dj^) . (93) 

It follows immediately that the conjugate momentum to the scalar field can 
be expressed as 

n{r],x) ^ —^'{r],x). (94) 

c 

It is convenient to Fourier expand the field ^{rj, x) over the basis of plane waves 
(therefore, here, we use explicitly the fact that the spacelike hypersurfaces are 
flat). This gives 

<^(.,.) = ^(^/dfe..(.)e-. (95) 

We have chosen to define the Fourier component with a factor l/a{ri) for 
future convenience. Since the scalar field is real, this last relation allows us to 
write /i^ = /i-fe. The next step consists in inserting the expression of <P{ri, x) 
into the action. This gives 



S = — [ dr] [ d^k 
2c J '7fl3+ 



Mfe Mfe + MfcMfc - 2— (/ifeA^fc + Mfc A*fc) 



a 



a'2 



---k ] ifik^^l + tJ*ktik) 



(96) 



Notice that the integral over the wavenumbers is calculated in half of the 
space in order to sum over independent variables only. Equipped with the 
Lagrangian in the momentum space (that, in the following, we denote by £), 
we can now go to the Haniiltonian formalism. The conjugate momentum to 
^fc is defined by the formula 

S^ 1 / , a' \ 

P'' = TT7 ^ -{fJ-k fJ-k ■ (97) 

Sn* c V a J 
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One can check that the definitions of the conjugate momenta in the real 
and Fourier spaces are consistent in the sense that they are linked by the 
(expected) expression 



n{r^,x)^j^^J dkpke^'^-^. (98) 

We see that the definition of the conjugate momentum pk as the derivative 
of the Lagrangian in the Fourier space with respect to /ij,* and not to /i^ is 
consistent with the expression of the momentum in the real space. Otherwise 
the momentum iT(r/,a;) in real space would have been expressed in terms of 
p^ instead oi p^. 

One can also check that the Lagrangian leads to the correct equation of 
motion. Since we have SC/Sfil = l/{2c)[-2n^.'y, + 2{n'^ - k'^)ti'^], the Eulcr- 
Lagrange equation d[SC/S^'^*]/di] — 6£./S^l. ~ reproduces the correct equa- 
tion of motion for the variable /i^, namely 

^^ ■ ''^' 'Mfe=0, (99) 



dyy^ \ a 

which is indeed the well-known result. 

We are now in a position where we can go to the Hamiltonian formalism. 
The Hamiltonian density, 7i, is defined by 

H = pk^i'k*+P*kfi'k-£, (100) 

and we obtain 

/ P \ a' 

'H = c[ pkpl + —ti.ktJ*k + — (Pfc/4 + PfeMfc) • (101) 

One can check that the Hamilton equations 

dr] opk a drj Ofik a c 

lead to the correct equation of motion given by Eq. (|99|l . 

As a preparation to canonical quantization, we now introduce the normal 
variable ak |19| defined by 

ak iv) = N{k)fik + ic—^pk , (103) 

k 

where, for the moment, the functions N{k) and M{k) are free but will be 
specified later on. In terms of the normal variables, the scalar field and its 
conjugate momentum can be expressed as 

1 1 r dk 

^(^' ^) - ^ (2^ y 2ivM [-('^^^"-^ + "^('^)^""-^] ' ^'''^ 

niv,^) -^.j^^^) [«.(.)e---«U.)e--] . (105) 
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We are now ready to quantize the system. So far, we were dealing with a 
relativistic field theory and only the constant c appeared in the equations. 
Now, the constant h, which fixes the amplitude of the fiuctuations shows 
up. Concretely, the quantization is carried out by requiring that <?(?y, x) and 
77(77, ^) become quantum operators satisfying the usual commutation relation, 
namely 

[S{r,,x),n{r],y)]^ihS^{x-y). (106) 

The normal variable akiji) is promoted to an operator Ckiji). We choose the 
commutation relation to be [cfc(r7), cj,(r/)] — CS{k ~p). In the last expression, 
C is a free dimensionless constant. Notice that the commutation relation is 
time- independent. Then, the expressions of N{k) and M(fc) are fully deter- 
mined. Let us see in more details how the calculation proceeds. The commu- 
tator is given by 



[<^(?7,a;),77(ry,y)] 



iC 



4c(27r)3 



d-^fe- 



k 



N{k)M{k) 



.Jk(x-y) 



-ik(x—y) 



(107) 



We see that, in order to produce a Dirac function S^{x~y) which is necessary 
in order to reproduce the relation given by Eq. (|106|) by integration of the 
exponentials, the term k/{NM) must be fc-independent. The link between 
the functions N{k) and AI{k) is therefore determined. Let us call D the term 
k/{NM). Then the result reads 



[<?(ry, x),n{r^, ?/)] - ^ x 2DS''{x - y) . 
4c 



(108) 



As a consequence we have CD = 2hc. As expected, the normalization is given 
by a combination of h and c. In the following, we will adopt the convenient 
choice C = 1. 

Everything has been fixed but the function N(k). This function is chosen 
by means of the following considerations. The energy of a scalar field is given 
by the formula 

E^ I Sx^/-^^)gp^ f d^xy/^^^ (S'^ + S'^^^S^JS^ (109) 



d'xv/^raj- 



S'^dr^dj^ 



110) 



where we have used the expression of the conjugate momentum. In this ex- 
pression, the determinant of the metric is the determinant of the spatial part 
of the metric (including the factor a). We can now insert the expression of the 
operators <!> and 77 in the above equation giving E. One finds 



E 



1 

2^ 



1 



d-^fc- 



A72(fc) 



(- 



CkC-k 



Ckcl 



J. 



J J 



+ 



N^k) 



[ckC-k + Ckcl + clck + Cfcctfc) 



(111) 
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Our criterion is to put "half of a quanta in each mode" . Technically, this 
means that we would like the energy to take the following suggestive form 



E^ J d^fe^^ [ckcl + clck) , (112) 

where w(r/) = kc/a{ri) is the physical frequency. We see that the only way to 
cancel the unnecessary terms in Eq. (|111|1 is to have N{k) = M(k). Together 
with the relation established previously, D — 2'hc, this gives N'^{k) = k/{2hc). 
As a consequence. The scalar field operator now reads 

Everything is now fixed. The expression of the scalar field operator contains 
no unspecified factor. Even the amplitude is fixed and is given by the factor 

We can now calculate the Hamiltonian operator. Using Eq. (|1()1|) one ob- 
tains 



H = i / d^fc 



hk (ckCk'' + c_fc^c_fc) - ih— (ckC-k + c-fe^cfc^) 
^ 'a 



, (114) 



where it is important to notice that the integral is calculated in R^ and not in 
R^^ . Let us analyze this Hamiltonian. The first term is the standard one and 
represents a collection of harmonic oscillators. The most interesting part is the 
second term. This term is responsible for the quantum creation of particles in 
curved spacetime. It can be viewed as an interacting term between the scalar 
field and the classical background. The coupling function ia' /a is proportional 
to the derivative of the scale factor and therefore vanishes in flat spacetime. 
From the structure of the interacting term, i.e. in particular the product of 
two creation operators for the mode k and — fc, we can also see that we have 
creation of pairs of quanta with opposite momenta during the cosmological 
expansion. 

We can now calculate the time evolution of the quantum operators (we 
are here in the Heisenberg picture). Everything is known if we can determine 
what the temporal behavior of the creation and annihilation behavior is. The 
temporal behavior is given by the Heisenberg equations which read 

-J— = -T[cfc,H], — — = --[cfeT,H]. (115) 

arj Ti arj n 

Inserting the expression of the Hamiltonian derived above, we arrive at the 
equations 

^=K+-cV, ^ = -K^--l.- (116) 

arj a ar] a 
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This system of equations can be solved by means of a Bogoliubov transfor- 
mation and the solution can be written as 



Ckiv) == Uk{r])ck{riini) + Vk{T])C-k'' iVini) : 
CkHv) ^ ul{v)ckHVini) + vl{f])C-k{mni) : 



(117) 

(118) 



where rji^i is a given initial time and where the functions Uk{i]) and Vkii]) 
satisfy the equations 



i — = kuk{ri) 



^—^kiv): i — -, = kvkir]) + i—Uk{r]) . 

a dr] a 



(119) 



In addition, these two functions must satisfy \uk\ 



Vk 



1 such that 



the commutation relation between the creation and annihilation operators is 
preserved in time. A very important property is the initial values of the two 
functions are fixed and, from the Bogoliubov transformation, read 



Ufc(?7ini) = 1, Vk{'nini)=0- 



(120) 



At this point, the next move is to establish the link between the formalism 
exposed above and the classical picture. For this purpose, it is interesting to 
establish the equation of motion obeyed by the function Uk + v^.. Straightfor- 
ward manipulations from Eqs. I|119l) lead to 



{Uk 



a 
a 



(uk+vD^O. 



(121) 



Therefore, the function Uk + v^ obeys the same equation as the variable ^k- 
This is to be expected since, using the Bogoliubov transformation, the scalar 
field operator can be re-written as 



<^(ry,x) 



he 



1 



dfc 



a(77) (2^)3/2 7 ^2fc 

+ (ul + Vk) iv)climni)e-'''-'' 



{uk +vl) (ry)cfc(r/ini)e 



ikx 



(122) 



If we are given a scale factor, we can now calculate completely the time evolu- 
tion of the perturbations by means of the formalism presented above. Let us 
stress again that the quantization procedure has completely fixed the overall 
amplitude of the field. Indeed, the field is normalized to vhc while the "mode 
function" Uk + f ^ has initially an amplitude of one. 

Let us now calculate the two-point correlation function in the vacuum 
state. One gets 



^{t], x)<P{ri,x + r) 



he 

4^ 



dfc sin kr 
k kr 



Uk+Vl 



a{m) 



(123) 
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If we assume that the scale factor is given by a power-law of the confornial 
time, a{ri) — ioi^v)^^^ j where (3 < —2 is a free a parameter and £q a constant 
with the dimension of a length, then the solution of Eq. (|121|) with the initial 
conditions given by Eqs. (|120|l reads 



{uk+vl){r]) = 



,j(fer,i„i-7r/3/2)^ 



v-fc^^ij_i/2 (-M > 



(124) 



where iJ ^^■' is a Hankel function of first kind. From, this solution, it is easy to 
calculate the spectrum on large angular scales {krj -^ 0) 



he . 



Uk + vl 



a{v) 



Tt-C /(/?), 4+2/3 



47r2 ^2 



(125) 



where /(/3) = tt~^ [2-'^-Pr {-(3 - 1/2)]^ In particular, if /3 = -2, this case 
corresponding to de Sitter spacetime for which the Hubble constant is strictly 
constant, one has -^o = c/ffj^f and the spectrum reads 



c V 27r 



(126) 



i.e. is scale-invariant (which means that it does not depend on the wavenum- 
ber). Of course, if /? 7^ —2 then the spectrum is scale dependent. 

The above result leads us to a first attempt to quantize cosmological per- 
turbations 120]. Using Eq. H73|l and taking into account the fact that uj^^^ ~ — 1 
and TL/ip' ~ K/{2e) (recall that e is the first slow-roll parameter) , one obtains 
for the spectrum of the "tracer" Cfc 



^C 



k\i^i-^k^ 



'5^ 



Sip 



(gi) 
k 



(127) 



Historically, the idea 



The question is now how should we calculate 

was to consider that the matter fluctuations (i.e. fluctuations in the scalar 
field) are quantized while the fiuctuations in the gravitational field remain 
classical. Based on this guess, one can used the trick which consists in replacing 



Sip 



(gi) 



Sep 



(gi) 



(128) 



or, in the Fourier space, 
of density perturbations 



Sep 



(gi) 



3/-2 



P. EE kX 



fiH'^^J {Att^ c) . This gives for the spectrum 



(129) 



As expected the three fundamental constants, G, c and h participate to the 
final expression. We have kept them in order to be able to trace back their 
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origin. The combination which appears here is the Planck time squared as 
it has to be since (k is a dimensionless quantity. In natural units, the above 
spectrum is just H^ /("ntm? ). Several remarks are in order at this point. 
Firstly, as we will see, this trick provides us with the exact result. Secondly, 
it seems is that there is no way to rigorously justify the replacement (|128ll . 
The reason is that matter, i.e. iST^j/, is treated quantum-mechanically, while 
geometry, i.e. 5G^u, is still considered to be classical, despite the fact that 



both are linked by the perturbed Einstein equations, ^G^jy 
could think that a semi-classical equation 



nST^iy. One 



6G 



fiiy 



5T^, 



(130) 



could do the job but in fact one easily realizes that this cannot be the case 
because ST^j^^ being linear in 5ip (it is of course quadratic in the scalar fields, 

o\ = 



but at linear order we have terms like tp'Sip'), we have in fact ( 5T^ 

due to (0 \5(p\ 0) = 0. Therefore, Eq. (|130|) is in fact inconsistent in the present 
context. Thirdly, it would be dangerous to base the physical interpretation of 
Eq. p29() on the above method arguing that it gives the correct result. Here, 
we emphasize that a convincing physical interpretation should be based on a 
consistent framework. We can try the following analogy. The correct equation 
for the energy levels of an Hydrogen atom. En oc 1/n^, has been obtained for 
the first time by means of the so-called Bohr's model. This model was devel- 
oped before a consistent framework for Quantum Mechanics become available. 
But, it is clear that, today, nobody would try to use Bohr's framework to in- 
terpret the formula for En- We are of the opinion that the situation for the 
cosmological perturbations is similar. Fourthly, the correct way to proceed is 
to treat the fluctuations in the geometry and in the scalar field on an equal 
footing. This amounts to "quantize" both sides of the Einstein equations and 
to write [H] 



5G, 



k6T,, 



(131) 



The consequence is of course of utmost importance: the metric operator h^^, 
should now be considered as a quantum operator, /i^t^ —^ h^^. In other words, 
we have now to deal with the quantum-mechanical nature of the gravitational 
field, i.e. with quantum gravity (at the linearized level). We now turn to this 
question. 



4.2 Quantization of Density Perturbations 

The total action of the system is given by 



IGttG 



d^a 



^B- 



1 



d^'a 



g'^'d^^d,^ + V{^) 



(132) 
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If we perturb this action up to second order in the metric perturbations and 
in the scalar field fluctuations (this is necessary if we want the variation of 
this action to reproduce the first order equations of motion) one finds, despite 
a very long and tedious calculation, that the result is delightfully simple, 
namely |S] 

1 " 



>dS 



2c 



d*x 



(«') 



l\2 



S'^^divdjV 



with 



v{ri, x) = a 



TL 



(133) 



(134) 



This is nothing but the action for a scalar field with a time-dependent mass. 
The constant G does not appear explicitly in the above action because it 
has been absorbed via the background Einstein equations. It is interesting to 
notice that the natural variable is not ^ neither C, but w. It is not a surprise that 
the system is characterized by a single quantity since gravitational fluctuations 
are described by <? and matter fluctuations by (5(^^8i) b^t are linked by the 
perturbed Einstein equations. Therefore, only one degree of freedom is left. 
The link between v and the "tracer" C, is given by 



c 



(135) 



Finally, the quantity Zg is given by Zg — y/nplaip' /TL = a^/e because, from 
the background Einstein equations, one has ni^p'Y — 27i^e. 

At this point, the procedure of quantization follows exactly the one pre- 
sented in the last subsection. The quantity v{rj,x) becomes a quantum oper- 



ator v{r],x) 
as 



<5(p(gO + (ip'/n) S 



, the expression of which can be written 



v{ri,x) 



he 



(27r)3/2 



d^fc 



(uk + vl){i])cki'nini)e 



ikx 



+ (.ul + Vk){v)<iiVinOe 



— ik-: 



(136) 



As announced, the fluctuations of the metric tensor are now quantized: techni- 
cally, the Bardeen potential ^{rj, k) is now a quantum operator ${ri, x) which 
explicitly appears into the expression of v{ri,x). Notice also that the un- 
certainty principle has completely fixed the overall amplitude of the quan- 
tum perturbations since the initial conditions are fixed by ^4,(771,11) = 1 and 
^fe(%ii) = 0. The equation of motion reads 



(uk+vl)" + 



Zc 



(ufc + vl)^0. 



(137) 



The physical meaning of the initial conditions are as follows: initially, we 
choose the state which is empty of "particles" from the point of view of a 
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local coiTLOving observer at the initial time r/inj. This state |0) is defined by 
Cfe|0) — 0. Since, due to the time dependence of the background, there is a 
nontrivial mixing between positive and negative frequencies, this state is in 
general not the vacuum at later times 

We are now in a position to calculate the power spectrum of the quantum 
operator (. One gets 



*/ N >/ NX S.CK /"*'°° dfc sinfcr ,,, ^,„ 

s 



Stt^z^ ./o ^ kr 



from which we easily deduce the expression of the power spectrum 

(139) 



hcK , 



87r2 



Uk + vl 



Zsiv) 



In order to compare this result with Eq. 11291) , we can evaluate the spectrum 
for power-law inflation where an exact solution of the equation of motion is 
available. Indeed, for a{ri) = tai—vY^^i the function e is a constant, hence 
one has z'^/ z^ — a" /a. This means that the Hankel function of Eq. (|124|l is 
also solution of Eq. H137() . Then, straightforward calculations show that 

If j3 is close to —2 then £o — c/H^^^ and one recovers exactly the result of 
Eq. (|129|l . As expected, the Planck length "naturally" appears in the above 
result. 



4.3 Quantization of Gravitational Waves 

The quantization of gravitational waves proceeds exactly along the same lines 
as before. Therefore, in this subsection, we only review briefly the main results. 
The starting point is the Einstein-Hilbert action that we expand to the second 
order. One gets (in natural units) |S] 

^''>SS =-^j [{h',y{hh)' - dk{h',)d''{hh)] a'{7^)d\. (141) 

In fact, the action can be re- written as 

^'^SS2^-j^ J2 U^x-g^'^d^h^d,h\ (142) 



where the quantity h'^{'q, x) is defined by 



2 

(^,-)-^(^ E /dfcA^;(^,fe)e*=^. (143) 
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Therefore, the action of gravitational waves is equivalent to the action of two 
decoupled scalar fields (corresponding to the two states of polarization) . One 
can then follow the method presented before. The quantum perturbed metric 
operator can be written as 



hij{r],x) 



40F 



1 



E 



dfc 



m,Av) (27r)3/2 J^^^ 



+{<* + vl)ivHHmnOe-"-' 



vUk) 



K + vinivHimmy' 



(144) 



where the function (u^ + w^ *)(??) obeys Eq. (|76(l . Finally, the two-point cor- 
relation function of the perturbed metric operator can be expressed as 



hij(rj,x)h'-^{ri,x + r) 



16 



Trmiyi'i^) Jo 



dfc sin fcr 
k kr 



k'\ul + vi*\\ 



(145) 



from which wc deduce that the power spectrum of the gravitational waves is 
given by 



fc3P;,(fc,77) = 



16 



k^ 



a{ri) 



(146) 



If we had decided to keep the standard units, the factor l/rn^ in the above 
result would have obviously read hG/c^, i.e. the Planck length squared. 



4.4 The Power Spectra in the Slow-roll Approximation 

We have established the expression of the scalar and tensor power spectra 
and calculated these quantities for power-law inflation. However, as discussed 
at the beginning of this review article, the most interesting physical situation 
occurs when the slow-roll approximation is valid. As discussed previously, the 
only thing we need to do in order to compute the spectrum is to solve the 
equation of a parametric oscillator, 



^,"+[k'-u{v)]^, = ^^" + 



iv) 



M = 0, 



(147) 



where ^ is Uk + vl either for scalar or tensor perturbations and the effective 
potential z'^/z^ or a"/a, i.e. z = z^ or z = a{r]). As already mentioned, 
on subhorizon or superhorizon scales, this equation can be solved regardless 
of the detailed form of the scale factor. The solutions are exp (— ifcry) and 
z{t]) + z{ri) J drz~^(r) respectively. However, in order to obtain a reliable 
solution, one also needs to know the form of the solution in the regime fc^ ~ 
U{r]), that is to say when the corresponding scales crossed out the horizon 
during inflation. 
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Let N^ (A) be the number of e-folds before the end of inflation at which 
the scale A exits the horizon. We have 



iV,(A) 



'■■'r 



log 



10 



H. 



-logio(5^ 



29 



xlnlO. (148) 



If we take the fiducial values K^, ~ lOi-^GeV and T^^ ~ A/,„j ~ lOi^-^GeV 
then A'"* ~ 60 for the Hubble scale today, i.e. A ~ £„, see also Fig. ^ A scale 
characterized by its wave-number k corresponds to an angle 9 on the celestial 
sphere of about k ~ l/{2£-^6). Given the present CMBR experiments, this 
means that we probe in fact the scales ^jj < A < 10~'^£jj. The smallest scale in 
this interval crossed out the horizon ~ 46 e-folds before the end of inflation. 
This means that the time taken by the the scales of astrophysical interest 
today to cross the horizon during inflation corresponds to AN ~ 7. Therefore, 
we need an accurate description of the effective potential U{ri) only during 7 
e-folds, see Fig. Q] 

In the slow-roll approximation, the effective potentials for scalar and tensor 
read at linear order 



Moreover, the equations of motion for e and 5 can be written as: 



de 



de 

AN 



2e{e~5) 



A5 



d(5 



2e(e-<5)-e. 



(149) 



(150) 



From these equations, one sees that typically O (e^) AN <C 1 for AN ~ 7 and, 
therefore, the slow-roll parameters can be considered as constant during the 
exit of the physical modes. This simplifies the problem drastically since then 
the equations of motion in the regime k^ ~ U{ri) can be solved in terms of 
Bessel functions whose orders depend on the slow-roll parameters. A detailed 
calculation can be found in Refs. j21| and, here, we just give the result 



k-^Pr 



k^Ph. 



i/2 



PI 

16H^ 



mn^ 



1 - 2 (C + 1) e - 2C (e - (5) - 2 (2e - ^) In 



k 



rC^: 



l-2(C + l)e-2eln — 



(151) 
(152) 



where C is a numerical constant, C ~ —0.73 and fc, a scale called the "pivot 
scale" . We see that the amplitude of the scalar power spectrum is given by 
a scale-invariant piece, H^ /{nem? ) that we had already guessed before, plus 
logarithmic corrections the amplitude of which is controlled by the slow-roll 
parameters, i.e. by the microphysics of inflation. It is important to notice 
that H is the value of the Hubble parameter during the 7 e-folds where the 
scales of astrophysical interest crossed out the horizon, see Fig. Q] As already 
mentioned at the end of Sec. (2.4) this can be different from the value of the 
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Hubble parameter at the beginning of inflation. The above remarks are also 
valid for tensor perturbations. The ratio of tensor over scalar is just given by 

This means that the gravitational are always sub-dominant and that, when 
we measure the CMBR anisotropies, we essentially see the scalar modes. This 
is rather unfortunate because this implies that one cannot measure the energy 
scale of inflation since the amplitude of the scalar power spectrum also depends 
on the slow-roll parameter e. Only an independent measure of the gravitational 
waves contribution could allow us to break this degeneracy. On the other hand, 
the spectral indexes are given by 



n, - 



dlnfc 



1 - 4e -f 2(5 , n„ 



k=k, din A: 



==-2e. (154) 



As expected, the power spectra are always close to scale invariance and the 
deviation from it is controlled by the magnitude of the two slow-roll parame- 
ters. Finally, at the next-to-leading order there is no running of the spectral 
indexes since they are in fact second order in the slow-roll parameters. 



5 Comparison with Observations 

In this section, we briefly discuss the impact of the recent Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) measurements on inflation [32]. We 
have seen previously that the presence of cosmological perturbations causes 
anisotropies in the CMBR (the Sachs- Wolfe effect) and we have established 
the link between 5T/T and the metric fluctuations, see Eqs. ifS^ and (|OT|l . 
The fact that the metric fluctuations are described by a quantum operator 
has an immediate consequence: 5T /T should be considered as a quantum op- 
erator as well. It is convenient to expand this operator on the celestial sphere, 
i.e. on the basis of spherical harmonics 

t~^ -t-oo ra^H. 

01 



T 



e)=^ ^ a^™y«,„(6',.^). (155) 



The next step is to calculate the two-point correlation function of temperature 
fluctuations. One gets 







^(61)^(62) 



+00 



e=2 



(2^_+l 
An 



0)^E^^^^''^^(cos^)' (156) 



where Pi is a Legendre polynomial and 7 is the angle between the two vectors 
ei and 62- The Cg 's are the multipole moments and have been measured with 
great accuracy by the WMAP experiment |22) . 
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A remark in passing is in order at this point. As a matter of fact, what 
has been measured by the WMAP sateUite is the correlation function 

f (e.)f (e.)) , (157) 

where the bracket denotes spatial average over the celestial sphere and not 
ensemble average as in Eq. H156|) . Going from one to another is not trivial and, 
in fact, involves profound questions which can even go as further as problems 
linked to the interpretation of Quantum Mechanics! (another related question 
is the problem of the "classicalization" of the quantum perturbations, see 
Refs. US]). In order to check that the predictions of Eq. H156|l are verified or 
not, one should repeat the measurement of the CMBR map many times and 
see whether the result converges toward the theoretical prediction. However, 
one cannot do that because we only have at our disposal one realization, i.e. 
one Universe or one CMBR map. Facing this situation, the usual strategy is to 
construct an unbiased estimator of the quantity that we want to measure (the 
correlation function or the multipole moments) with the minimum possible 
variance so that it is very probable that the outcome of one realization is 
closed to the mean value [23. Unfortunately, the variance cannot be zero (in 
this case only one realization would be enough to estimate the result) and one 
can show that this is linked to the fact that a stochastic process on a sphere 
cannot be ergodic |21]. This variance is called the "cosmic variance" and is 
generally large on large scales. More details on this question can be found for 
instance in Ref. j24| . 

On large scales, i.e. for small ^, one can use Eq. 1)91(1 to find an explicit 
expression of the multipole moments. One gets 

A f+OC 17 

^'^2EJ„ T^IW^'^C, ^«20, (158) 

where ji is a spherical Bessel function of order £. Using Eq. ((151|l for density 
perturbations (since they are dominant) and neglecting the logarithmic cor- 
rections (which amounts to consider that the spectrum is scale-invariant), we 

obtain 

2if2 1 

25em2^ £{£ + 1) 

Therefore, a scale invariant spectrum implies that, on large scales, the quan- 
tity £(£ + l)Ci is a constant. In order to calculate the inflationary multipole 
moments C( for any £ one must use a numerical code, for instance the CAMB 
code 25 . Typically, one gets a plateau and then acoustic oscillations. Here 
we do not treat this question but the details can be found in Ref. [TH| . 

The sate llites COBE and WMAP have measured the quantity Q/T = 
a/5CV(470 where T ~ 2.7K and have found Q ~ 18 x IQ-^K. Moreover, 
recent analysis |2S| of the WMAP data have been able to put a constraint 
on the value of the slow-roll parameter e. It was found that e < 0.032. This 



Q = 7r7:-^77^^T' ^<20. (159) 
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allows us to put a constraint on the Hubble parameter at horizon crossing. 

One finds 

i?2 n2 



= 607re^ 






< 1.6 X 10" 



(160) 



This also puts a constraint on the amount of gravitational waves. In Ref. |26j . 
the following result has been obtained 






<0.3. 



(161) 



that is to say the contribution of gravitational waves is already constrained 
to be less than 30% of the total contribution. 

We conclude this part by a summary of the main observational predictions 
of single field inflation: (i) The universe is spatially flat: X?o = 1 ± 10~^; (ii) 
The spectrum of density perturbations is scale invariant (Harrison-Zeldovich 
spectrum) plus logarithmic corrections which are model dependent, i.e. n^ = 
1 + 0{e, S); (iii) There is a nearly scale invariant background of gravitational 
waves, i.e. n^ = 0(e); (iv) The statistical properties of the CMB anisotropics 
are Gaussian, i.e. everything is characterized by the power spectrum and we 
have the following properties 





3 ( 




0, etc 



(162) 

This conclusion comes from the fact that the quantum state of the perturba- 
tions is the vacuum, the "wave function" of which is a Gaussian; (v) Gravi- 
tational waves are sub-dominant and there exists a consistency check relating 
the importance of gravitational waves with respect to scalar density on one 
hand to the tensor spectral index on the other hand. This relation reads 



■pi" - "/2(^, ^cdm, ^A, ■ ■ O^T : 



(163) 



where the function /2 is /2 — 5 for the concordance model (i.e. the cold dark 
matter model plus dark energy which seems to fit best the data at the time 
of writing); (vi) There are oscillations in the power spectrum. Although this 
conclusion is also based on the physics of the transfer function, the fact that 
the perturbations are generated in a coherent manner plays a crucial role for 
the survival of the acoustic peaks, see Ref. [77) . 



6 The Trans- Planckian Problem of Inflation 



We have seen that the CMBR anisotropics are, if the infiation theory turns 
out to be correct, an observable signature of quantum gravity. However, as it is 
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clear from the previous considerations, the CMBR anisotropics originate from 
a regime where the quantization of the gravitational field is carried out in the 
standard manner. In fact, the situation is similar to the Hawking radiation. In 
this last case, we have a quantum field living in a classical background. In the 
present context, we also have a field h^^ijj^x) living in the classical FLRW 
Universe. Of course, the main difference is that, in the case of inflation, the 
quantized test field is the perturbed metric, i.e. is the gravitational field itself 
(at least the small excitations of the gravitational field around a classical back- 
ground) contrary to the Hawking effect where the field is just a scalar field: 
this is why, conceptually, the Hawking effect does not involve quantum grav- 
ity while the theory of cosmological perturbations does. Nevertheless, from 
the pure technical point of view, we have just used the techniques of ordinary 
quantum field theory in curved space-time. In this section, we suggest that the 
CMBR anisotropics could also carry some signatures of quantum gravity but, 
this time, originating from the non perturbative regime |2H1- Obviously, the 
price to pay is that the following considerations are much more speculative 
than the rest of this review article but the hope is to learn about quantum 
gravity, maybe in the non-linear regime. Therefore, it seems that the potential 
reward is worth the speculation. 

The inflationary trans-Planckian issue is based on a very simple re- 
mark |28j . If we assume a model, for instance a potential of the type given by 
Eq. (|2I|I (here, we choose n = 4 to be concrete), then one can calculate the 
coupling constant A„. For this purpose, it is convenient to express everything 
in terms of iV*, the number of e- folds before the end of inflation at which the 
modes crossed out the Hubble radius, see Eq. (|148|) . The corresponding value 
of the inflaton field is given by ifl = m'^^{N^, + 1)/tt. Therefore, the Hubble 
parameter can be expressed as H^ = X4m'^^{ip,,/mpJ'^ = \4m'^^{N^, + l)^/7r^. 
Finally, since the slow-roll parameter e is given by e = {N^, + 1)^^, one arrives 
at 

J^ = ^X4N, + lf. (164) 

em'' TT^ 
PI 

The scale of inflation only enters the above equation through N^ and the 
corresponding dependence is logarithmic, see Eq. H148|l hence very mild. One 
can thus use this formula to determine the coupling constant almost indepen- 
dently of H.^^.. Using Eq. H159|l for i = 2 and the link between Q and C2, one 
finds that A4 ~ 10~^^, where we have used A^* ~ 60. As already mentioned, 
this means that the total number of e- folds is huge, TV^, ~ 4.9 x 10^. As a re- 
sult, the Hubble radius today, £-^ — 10^^i?p, {h = 0.5), where ^pj is the Planck 
length, was equal to ~ e~^° £p[ ~ xo~'*-^^^° £-p^ at the beginning of inflation, 
i.e., very well below the Planck length! 

One can view the problem differently and ask how many e-folds before the 
end of inflation a given scale was equal to the Planck length. The answer can 
be easily calculated from Eq. (|148|l and reads 
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TVp, (A) = N, (A) - logio f:^) X In 10 (165) 
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29 - logio ( ^^ 



m^ 



xlnlO. (166) 



If one takes the fiducial values H.^^ ~ IQi'^GeV, T^^ = M,„f ~ IQi^-^GeV, one 
finds that the Planckian region was reached only 11 e- folds before the modes 
crossed out the horizon during inflation, see Fig.^ For instance, for the mode 
A = fjj, this means 70 e-folds before the end of inflation. Of course, if the 
scale of inflation is smaller, then the number of e-folds before the exit of the 
Planckian region and the exit of the horizon can be bigger. 

The following point should also be emphasized. At the time at which the 
modes of astrophysical interest today exit the Planckian region, the value 
of the Hubble parameter is generically well-below the Planckian mass. This 
means that the use of a classical FLRW background is well justified. The trans- 
Planckian problem concerns only the fluctuations and has to do with the fine 
structure of the Universe or with the "Planckian foam" but does necessitate a 
full quantum gravity description of the evolution of the underlying manifold 
(for instance, one does not need quantum cosmology). 

Having in mind the above considerations, the trans-Planckian problem 
of inflation consists in the following (2£j. It is likely that the framework of 
standard quantum field theory described in the previous section and used in 
order to establish what the predictions of inflation are breaks down when 
the modes under consideration have a wavelength smaller than the Planck 
length. Therefore, there is the danger that the so far successful predictions 
of inflation are in fact based on a theory used outside its domain of validity. 
In other words, there is the problem that the predictions of inflation could 
in fact depend on physics on scales shorter than the Planck length, a physics 
which is clearly largely unknown. 

Is it really so? In trying to answer this question we immediately face the 
problem that the trans-Planckian physics is presently unknown and that, as 
a consequence, it is a priori impossible to study its influence on the infla- 
tionary predictions. To circumvent this difficulty, one studies the robustness 
of inflationary predictions to ad-hoc ( "reasonable" ) changes in the standard 
quantum field theory framework supposed to mimic the modifications caused 
by the actual theory of quantum gravity. If the predictions are robust to some 
reasonable changes, then there is the hope that they will be robust to the 
modifications induced by the true theory of quantum gravity. On the other 
hand, if the predictions are not robust, the knowledge of the exact theory 
seems to be required in order to predict exactly what the changes are. The 
next question is of course which kind of modifications can we introduce in the 
theory in order to test its robustness? Many proposals have been made and 
discussed recently in the literature j28l 1291 ^^ \'M\ . Here, we concentrate on 
two possibilities: the modified dispersion relation and the so-called "minimal" 
approach. 
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6.1 Modified Dispersion Relations 

Let us start with the modified dispersion relations. The term k^ in Eq. H147|) 
originates from the use of the standard dispersion relation uJ^^^y^ = fcphya ■ In 
condensed matter physics, it is known that the dispersion relation starts de- 
parting from the linear relation tu = k on scales of the order of the atomic 
separation: the mode feels the granular nature of matter. In the same way, one 
can expect the dispersion relation to change when the mode starts feeling the 
discreteness of space-time on scales of the order of the Planck (string) length. 
Therefore, our method is to replace the linear dispersion relation ui ^ = k , 

' ^ ^ phys phys 

by a non standard dispersion relation lu ^ = tj ^ (k), this non linear relation 
having of course the property that Wphy^ — ^phys fo^" k <^ k^ where k^ is a 
new scale introduced in the theory which could be, for instance the string 
scale. In the context of cosmology, this amounts to replacing the square of the 
comoving wavenumber fc^ with 



k' ^ kMk,v)^a^iv)tol 



a{ri) 



(167) 



Therefore, this implies that we now deal with a time-dependent dispersion 
relation, a result first obtained in Ref. [2H| ■ As a consequence, the equation of 
motion H147|l now takes the form 



M 



klsik,!!)- — 



M 



0. 



(168) 



The effect of the new physics is to change the time-dependent frequency iolk, rf) 
of the parametric oscillator. Let us remark that a more rigorous derivation of 
this equation, based on a variational principle, has been provided in Ref. |3U| . 
Then, the only question is whether the fact that we now have a new 
time-dependent frequency can modify the spectrum /c'^|/ip or not? As we now 
demonstrate, this depends on whether the evolution of the modes is adiabatic 
or not in the trans-Planckian region. Indeed, if the dynamics is adiabatic 
throughout (in particular if the z" / z term is negligible), the WKB approxi- 
mation holds and the solution is always given by 
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(169) 



where ryjni is some initial time. Therefore, if we start with a positive frequency 
solution only and uses this solution, one finds that no negative frequency 
solution appears. Deep in the region where fcoff — k, i.e. for k <ti k^, the 
solution becomes 



Kv) 



1 



exp 



—ikrj — i j kcs{k,T)dT 

r;i„i 



(170) 
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where rji is the time at which k^s ~ k. Up to an "accumulated" phase which 
will disappear when we calculate the modulus |/xp, we recover the standard 
vacuum solution e~*'^''/v2fc and hence the standard spectrum. We have thus 
identified the criterion which controls whether the spectrum will be changed 
or not: in order to get a modification, the dispersion relation in the trans- 
Planckian region must be such that the WKB approximation is violated. This 
constrains the shape of the modified dispersion relation. It is possible to give 
the conditions for violation of the WKB approximation. Given an equation of 
the form fi" +uj^ijl = (in the present context, one has w^ = fc^g — z" / z), the 
WKB approximation is valid if the following quantity is small in the trans- 
Planckian region [S^] 

^|«1, (171) 

where Q is defined by the following expression Q — 3(a;')^/(4a;^) — uj" /{2lo). 
Then, one can insert in the previous expression one's favorite dispersion rela- 
tion ans see whether this leads to a new spectrum. This has been done recently 
in the literature, see Refs. [2^1 ■ For instance, one can show that the dispersion 
relations introduced in Refs. \i'i\ I34j do not lead to any modification. An ex- 
ample where modifications are present has been studied in Ref. j30j . However, 
it remains to be studied whether this can be made compatible with other 
studies on the subject, in particular those using astrophysical observations 
to constraint the deviations from the law ui = k [2Sj. Rather than studying 
these examples in great details, we now turn to a new way of modeling the 
trans-Planckian regime. 

6.2 The Minimal Approach 

Modifying the dispersion relation is equivalent to changing the form of the 
equation of motion for the perturbations. The minimal approach consists in 
working with the same equation of motion (with a standard dispersion relation 
hence the name "minimal approach") but with modified initial conditions. 
For a given Fourier mode, the initial conditions are fixed when the mode 
emerges from the trans-Planckian region, i.e. when its wavelength becomes 
equal to a new fundamental characteristic scale i^-. = l/fc^. The time rjk of 
mode "appearance" with comoving wavenumber k, can be computed from the 
condition 

9-77 Ptt 

X{Vk) = j^aM=i,^—, (172) 

which implies that rjk is a function of k. This has to be compared with the 
standard inflationary calculations where the initial time is taken to be r]k — 
— oo for any Fourier mode k and where, in a certain sense, the initial time does 
not depend on k. Then, a crucial question is in which state the Fourier mode 
is created at the time rjk (here, we cannot take the limit ki] ~* ^oo anymore). 
The only requirement is that, if we send the new scale M^ to infinity (i.e. 
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if there is no trans-Planckian region), then one must recover the standard 
WKB vacuum. Therefore, the most general parametrization of these initial 
conditions read 



IJ-iVk) 



Cfc + dk 40r 



IJ-'iVk) 



±i 



^s,T('7fe) 4V7r(cfe -4) 



(173) 



where the coefficients Ck and dk are a priori two arbitrary complex numbers 
satisfying the condition |cfcp — \dk\'^ — 1 and which can be expanded as 
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(174) 



where ctq = H/M^. When M^ is sent to infinity then ctq ^ 0, Cfc = 1, dfc = 
and, indeed, we recover the standard vacuum. Since there are two energy 
scales in the problem, namely the Hubble parameter H during inflation and 
the new scale M^, it is natural that the final result is expressed in terms of 
their ratio H/M^ , which is typically a small parameter. The parameters x and 
y are considered as free parameters that are not fixed by any existing well- 
established theories except, as already mentioned above, that they should be 
such that the relation jc^p — jdfcp = 1 is satisfied. One easily shows that this 
implies y + y* = at leading order in ctq. Expanding everything in terms of 
ao, one arrives at ^T] 
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where 95 is the argument of the complex number x, i.e x = \x\e^'^ . These ex- 
pressions should be compared with Eqs. I|151l) and (|152|l . The effect of the 
trans-Planckian corrections is clear: superimposed oscillations in the power 
spectra have appeared. The magnitude of the trans-Planckian corrections are 
linear in the parameter a^ and their amplitude is given by |a;|cro- The wave- 
length of the oscillations can be expressed as Ak/k = <7o7r/e. 

The above calculation provides us with an explicit example where the 
observational predictions of inflation arc modified by the trans-Planckian 
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physics. Let us now study this question in more details. Using Eq. H91|l . one 
can evaluate the modifications of the multipoles moments caused by the trans- 
Planckian corrections. In the limit e/co ^ I., one gets [361 137| 



(176) 
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This expression should be compared with Eq. (|159|l . The oscillations in the 
power spectra are transfered to the multipole moments, at least at relatively 
small scales. At large i, or for not too small values of ctq, the above equation 
quickly becomes invalid and an accurate estimation can be made only with 
the help of numerical calculations. The result in plotted in Fig. O for the 
temperature fluctuations but also for the polarization, for details see Ref. [3^1 
I37j . In those references, a detailed comparison of the trans-Planckian signal 
with the recently released high accuracy WMAP data has been performed. 
The main result is that, with the oscillations taken into account, it is possible 
to decrease the x^ significantly. Instead of x^ — 1431 for 1342 degrees of 
freedom for the standard slow-roll power spectra, one now obtains x^ — 1420 
for 1340 degrees of freedom, i.e. A}^ ~ 10 compared to WMAP one. The 
reason for such an important improvement of the y^ is due to the presence 
of the oscillations which permit a better fit of the cosmic variance outliers 
at small scales. The main question is of course the statistical significance of 
this result. In Ref. [201137], the so-called F-test has been used and indicates 
that the result is significant. However, it is clear that other statistical tests, a 
complete exploration of the parameter space and, of course, new data, should 
be used before one can really conclude that superimposed oscillations are 
really present in the CMBR multipole moments. A fair description of the 
present situation is that there seems to be a hint for an interesting feature 
in the CMBR data and that, maybe, this feature is a signature of very high 
energy physics (it is clear that the oscillations, if their presence is confirmed, 
could have another physical origin). 

Finally, we would like to conclude by a comment on the back-reaction 
problem. This question is crucial for the consistency of the approach used 
before. It is clear that the energy density of the perturbations must be smaller 
or equal than that of the inflationary background. This leads to the condition 
\x\ < vSTTTOpi /A/^ which amounts to 

2 

\x\ao < 10^ X -^ . (177) 

It is important to emphasize that the above constraint is only a sufficient 
condition, but by no means, unless proved otherwise, a necessary condition. 
In general, this constraint is difficult to satisfy. Some of the best fits described 
above suffer from this back-reaction problem, see Ref. |3()l I3fil I37j . In fact, the 
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Fig. 3. Angular TT, TE and EE power spectra for two different trans-Planckian 
models, one with low frequency (LF) superimposed oscillations, the other with high 
frequency (HF) oscillations, for details see Refs. i36u37^ . A zoom of the temperature 
multipole moments in the first Doppler peak region is also shown (black curve) and 
compared with the standard slow-roll prediction (blue curve) calculated with the 
same cosmological parameters. 
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above formula expresses a generic difficulty of the trans-Planckian question, 
this difficulty being present regardless of the approach used in order to model 
the new physics. The presence of trans-Planckian corrections means the pres- 
ence of particles (with respect to the standard vacuum) the energy density 
of which is very easily of the order of the background energy density. On the 
other hand, if we try to satisfy the back-reaction constraint then the signal 
very easily becomes tiny and, hence, non observable. A major advance, which 
would allow us to escape the previous vicious circle, would be to calculate ex- 
plicitly the effect of the back-reaction. Unfortunately, for the moment, this is 
still an open question and more work is required to tackle this very important 
task. 
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